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In the past few years, coupling strengths between light and mechanical motion in optomechanical 
setups have improved by orders of magnitude. Here we show that, in the standard setup under con- 
tinuous laser illumination, the steady state of the mechanical oscillator can develop a non-classical, 
strongly negative Wigner density if the optomechanical coupling is large at the single-photon level. 
Because of its robustness, such a Wigner density can be mapped using optical homodyne tomogra- 
phy. These features are observed near the onset of the instability towards self-induced oscillations. 
We show that there are also distinct signatures in the photon-photon correlation function g^-^^ (t) in 
that regime, including oscillations decaying on a time scale not only much longer than the optical 
cavity decay time, but even longer than the mechanical decay time. 



By coupling optical and mechanical degrees of free- 
dom, the emerging field of optomechanics provides excit- 
ing new opportunities to study the quantum mechanical 
behaviour of macroscopic objects (for reviews see 
Recent optomechanical cooling experiments are success- 
fully bringing nanomechanical oscillators into their quan- 
tum mechanical ground state [31 13] . The same optome- 
chanical coupling also promises the possibility of back- 
action evading quadrature measurements of the result- 
ing mechanical quantum states with the help of the light 
field j4, J), • For a reproducible and persistent quantum 
state, such measurements would result in an experimen- 
tal determination of its full Wigner density via tomog- 
raphy, similar to what has been achieved in microscopic 
systems, for single ions or photons [2l[8]. The recent ad- 
vances in fabricating optomechanical devices have dras- 
tically improved coupling parameters, e.g. for optome- 
chanical crystals [S], in microwave setups [3], and other 
devices like GaAs disks fTU^ . It will likely be possible rela- 
tively soon to achieve optomechanical coupling strengths 
go at the single-photon level that are comparable to the 
optical cavity decay rate k, a feat that has already been 
achieved in cold atom optomechanical systems [TTJ [T^ . 
This regime of strongly nonlinear quantum optomechan- 
ics promise to pave the way towards generating and de- 
tecting novel quantum states in optomechanical systems. 
It is currently only beginning to be explored theoreti- 
cally P"5HT5] , although very early work already discussed 
quantum optomechanical effects in the (unrealistic) ab- 
sence of any dissipation [THl [IT] . 

In the classical regime, nonlinear dynamics is known to 
occur when the system is driven by a blue detuned laser. 



When the input laser power crosses a certain threshold, 
the mechanical oscillator will undergo a Hopf bifurcation 
and start self-induced mechanical oscillations, a mechan- 
ical analogue of lasing action [T5U21j . The quantum dy- 
namics of this regime has first been studied in [T3] , and 
there is interesting synchronization behaviour for arrays 
of coupled oscillators of this type [22 . 

In this paper, we show that, for strong optomechanical 
couplings, large laser driving and sufficiently low temper- 
atures, a non-classical state of the mechanical oscillator 
with strongly negative Wigner density can be produced 
around the onset of self-induced oscillations. Because the 
state is time- independent, one may use back-action evad- 
ing homodyne tomography to experimentally reconstruct 
its non-classical Wigner density. 

In addition, we propose to use the two-point photon 
correlation function g^^\t) as an experimentally conve- 
nient probe for the peculiar quantum dynamics near the 
bifurcation. We identify two distinct signatures that en- 
able experimentalists to reliably detect the onset and 
growth of the self-induced oscillation. We provide an 
explanation of the non-classical decay of g'^'^^ [t) in both 
the red and blue-detuned regime. 

A standard optomechanical system can be described 
by the following Hamiltonian: 

H = h{-A+goiP +b))a^ a+hujMb^b+haL{a'' +a)+Hdiss- 

(1) 

Here a/b are the operators for the photon/phonon modes, 
ujm is the mechanical frequency and is the laser driv- 
ing amplitude. A = lul — ujc is the detuning of the laser 
frequency from the resonant frequency of cavity at the 
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FIG. 1: The quantum self-induced oscillations in an op- 
tomechanical system. The laser power is held constant and 
the laser detuning A is being adjusted. The bottom right 
plot shows that, above a threshold, the phonon number ntof 
the mechanical oscillator rises abruptly. The other three 
plots show the evolution of the mechanical Wigner den- 
sity as the self-induced oscillation starts, with bars showing 
the reference sizes of ixzpp and 4pzpF- xzpf /pzpf are 
the zero-point fluctuation of the mechanical oscillator's posi- 
tion/momentum. Here go — 0.36, k — 0.3, F = 0.00147, Q_t = 
0.311, in units of um. The photon number in the cavity is 
Ua « 0.1—0.7 when -1< A < 0. 



mechanical displacement 0. go describes the strength of 
the optomechanical coupling at the single-photon level. 

When the dissipative terms in i?diss are taken into 
account, the density matrix p of the combined photon- 
phonon system evolves according to the quantum master 
equation: 



dp ^r., [H,p] 



TV[b,p\ + KD[a,p\. 



(2) 



Here C is the quantum Liouville operator describing the 
time evolution of the density matrix p, where we incorpo- 
rate dissipation in the photon/phonon subsystems with 
decay rates k and F, respectively. The standard Lindblad 
term is given by I?[d, p] ^ OpO^ ~\{d'^dp + pd'^d) [30]. 

In this paper we are interested in the steady state so- 
lution of Eq. [2j where all the transient dynamics has died 
out. This is obtained numerically by finding the density 
matrix satisfying C[p\ = 0. Due to this dynamical sta- 
bility this state is thus ideal for making homodyne mea- 
surements of its mechanical Wigner density, in contrast 
to transient scenarios. 

Specifically we are interested in the mechanical Wigner 
density Wm{x,p) = ;^ /^^.(x - ylpwl^: + yle^W^dy, 
where pM is the mechanical density matrix, obtained by 
tracing out the optical degree of freedom from p. The 
Wigner density is the quantum analog of the classical 
Liouville phase space probability density. A negative 
Wigner density is a strong signature of a non-classical 
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FIG. 2: Left Top: An example of a negative mechanical 
Wigner density, in this case at A = 0. For other parame- 
ters see Fig. [l] Left Bottom: The negative Wigner density 
state A = Qiblue) shows narrower phonon number distribu- 
tion than the near coherent state A = —0.25 (red). Right: 
Two parameter regions where there are significant areas of 
negative Wigner densities. The "quantum parameter" is de- 
fined as = So/k. We vary go and A, while keeping a, k and 
F the same as in Fig. [l] 



state. Early investigations [13] of the optomechanical in- 
stability in the regime around go ^ k had not turned up 
nonclassical states. 

In Fig[T] we show the overall properties of the steady 
state solutions. As we increase the laser detuning while 
keeping the input laser power fixed in points A B — )> C, 
the phonon number in the mechanical oscillator rises 
sharply, signaling the start of the quantum self-induced 
oscillation. This is also reflected in the mechanical 
Wigner density Wm{xtP). Before the self- induced oscilla- 
tion (point "A"), the oscillator stays in its ground state, 
with the Wigner density a simple Gaussian. The Wigner 
density starts to broaden just below the threshold, as the 
susceptibility of the system diverges and quantum fluc- 
tuations are strongly amplified (point "B"). The steady 
state density matrix of a periodically evolving system can 
be understood as the average of the Wigner density over 
one period, if the oscillation phase is random. In the 
regime of self-induced oscillations, a coherent state ex- 
ecutes circular motion in phase space while keeping its 
shape, which is the picture presented at point "C" . 

However, such a simple picture of a coherent state ex- 
ecuting a semi-classical motion is inadequate for an op- 
tomechanical system with ~ /t, i.e. when one ap- 
proaches the optomechanical instability deep in the quan- 
tum regime. In such a system, we observe that for a range 
of detuning A and laser driving Q!l, the mechanical self- 
induced oscillation produces strongly non-classical states 
with large negative areas in the Wigner density. This 
can be seen in the example of Fig. [2] left panel. Nega- 
tive rims develop at amplitudes slightly smaller than the 
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average amplitude of oscillation. The bottom left panel 
of Fig. [2] shows that the phonon number distribution of 
this state has a reduced variance (sub-Poissonian) , and 
is thus closer to a single Fock state as compared with a 
coherent state [Mj . 

The right panel of Fig. [2] maps out the regions in pa- 
rameter space where negative Wigner densities occur. 
This 'phase diagram' is shown as a function of the "quan- 
tum parameter" ^ (see [13] for more about that param- 
eter) and of the laser detuning A, at a fixed value of the 
laser driving strength. It has been obtained by solving 
for the steady state of the optomechanical system under 
constant illumination, and the Wigner density is consid- 
ered as nonclassical if a sufficiently large area turns out 
to be negative [3T] The numerical results shown here in- 
dicate that, starting at a "quantum parameter" ^ = 0.8, 
the negative Wigner density states appear around detun- 
ing A = 0, and a second negative Wigner density region 
opens up at ^ = 1.6, around A/cjjvf = 0.9 at the first 
blue sideband, where the instability is driven efficiently. 
The phonon number distribution displays a pronounced 
narrowing, getting closer to a single mechanical Fock 
state. However, we find that still many photon/phonon 
levels are involved in the dynamics in the regime consid- 
ered here, and there seems to be no simple explanation 
starting from a discrete energy level scheme. 

These steady-state non-classical Wigner densities could 
be reconstructed via optomechanical QND quadrature 
detection |S] and subsequent quantum state tomogra- 
phy |26j . This merely involves illumination with an- 
other amplitude-modulated laser beam for read-out, as 
explained in ^ . When observed, these would provide an 
accessible example of non-classical states in a fabricated 
mesoscopic mechanical object. To date, there has been 
no experimental observation of non-classical Wigner den- 
sities in the domain of micro- or nanomechanical struc- 
tures. The experiment that came closest to that goal, 
and in the process did produce nonclassical mechanical 
Fock states, employed a complex multi-layered supercon- 
ducting circuit with piezoelectric coupling to a supercon- 
ducting qubit and ultrafast pulse sequences [57] . Fur- 
thermore in their setup the resonator lifetime is too short 
to permit the reconstruction of the full Wigner density. 
By contrast, once optomechanical parameters can be im- 
proved to reach the single-photon strong coupling regime, 
the scheme discussed here would be relatively straight- 
forward, being based on continuous laser illumination of 
an optomechanical setup whose fabrication is much less 
complex as it involves only one material. And the non- 
classical states obtained are long-living and amicable to 
quantum homodyne Wigner density reconstruction. In 
addition, there is the possibility that the parameters re- 
quired here may be reached in cold atom optomechanical 
setups [TT][T2]. 

The full reconstruction of the mechanical Wigner den- 
sity of an optomechanical system in the nonlinear quan- 
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FIG. 3: Time-dependence of photon-photon correlations near 
the regime of quantum optomechanical oscillations. "A,B,C" 
labels the same parameters as in Fig.[l] Apart from observing 
photon bunching and oscillations of the correlation function 
at the mechanical period, these plots show that there is a re- 
markably slow long-term decay near the onset of self-induced 
oscillations, e.g. at point "B" (see main text). The insets 
show that the fully developed oscillation beyond the thresh- 
old, at point "C", leads to higher harmonics. 

tum regime is obviously an enticing and challenging goal. 
Nevertheless, it requires repeated measurements of many 
identically prepared quantum states. It will be helpful 
to have other means of optically probing the quantum 
dynamics of the system near the relevant regime around 
the onset of the instability. When the system itself is 
in a steady state, the expectation values of physical ob- 
servables do not change in time. A suitable probe for 
the dynamics is then provided by the photon two-point 
correlation function: 

5(2)(t) = i44+!^I±!^. (3) 
((aU.))2 

The (• • • ) denotes the statistical average over p in Eq. [2] 
In steady state, g^^^^ does not depend on the initial time 
T. Photon correlations are readily accessible in nowa- 
days quantum optics experiments with single-photon de- 
tectors (e.g. using a Hanbury-Brown Twiss setup), and 
they have been successfully employed to characterize the 
change of photonic statistics upon transmission through 
nonlinear systems. The most important example is pho- 
ton anti-bunching in the resonance fluorescence of single 
photon emitters, which has recently also been predicted 
to occur in optomechanical systems for sufficiently strong 
coupling [13]. 

As can be seen in Fig. [3] there are clear signatures in 
the photon correlations when the quantum self-induced 
oscillations start (around point "B"). 

In particular, g^'^\t) persists at some value above unity 
over a very long time (middle panel. Fig. [3]). It can 
be proven that as long as the steady state of the sys- 
tem is not degenerate, we always have g'^'^\t) — ?► 1 -I- 
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FIG. 4: Quantifying the slow approach of g''^^ (i) — >■ 1 near the 
onset of the quantum self-induced oscillations, as observed in 
Fig.[| 5*2) (i) _ 1 obeys 

an exponential decay e in the 
long-time limit t ^ oo. The inset shows the decay time Tg 
peaking toward very large values around point B, i.e. A = 
0.6cja/. 



aexp i~t/Tg) in the long-time limit i — >■ oo. Here the de- 
cay rate l/r^ — Re(Ai), where Ai is the eigenvalue of the 
Liouville operator C in Eq. [2] with the largest non-zero 
real part, characterizing the slowest decay in the system. 
This signature can be examined by plotting \n{g^'^^ {t) — l) 
to extract Tg, which indeed agrees with the Ai obtained 
numerically from £ (see Fig. |4| . As can be seen in the 
inset, Tg rises strongly around the start of the quantum 
self-induced oscillation (point B), resulting in a long-term 
persistence of g^^^^ (t) above 1 (although it finally decays 
to 1). This is connected to the fact that classically the 
overall mechanical damping rate goes to zero near the 
Hopf bifurcation [19j . 

The second signature in g^^' is the appearance of higher 
harmonics when the quantum self-induced oscillation is 
fully developed (see insets of Fig. [s]). Both the long- 
term decay of g^'^\t) and the higher harmonics appear- 
ing above the bifurcation threshold can best be under- 
stood in a semiclassical picture. To this end, we ap- 
proximate the photon correlator via the classical inten- 
sity correlator, (|Q;(i -I- r)p|a(T)p)^. The light amplitude 
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is modulated harmonically by 



the mechanical oscillations, as detailed in ]T5]. However, 
to obtain the decay of the resulting oscillations in the 
correlator, we have to take into account the mechanical 
phase diffusion induced by the radiation pressure shot 
noise. This was considered in pS], but here we discuss 
another, quite straightforward analytical approach: 

<50(t) = ^^^^jl cos(0(t')), (4) 

Var((50(t)) = ..^ ^ {Sff{i^m) + Sff{-ojm)) ■ 

(mwM^j 4 

Here Sff is the force noise spectrum (see [Mj). We 



FIG. 5: Semiclassical approximation for photon correlations 
in the mechanical lasing regime of optomechanics. By adding 
the mechanical phase-diffusion to the classical light field dy- 
namics (red), one can qualitatively account for the slow long- 
time decay of the photon correlator g'-^-' (t) in the full quantum 
simulation (blue) . The lowest panel compares the classical so- 
lution without phase diffusion (here A = ujm)- 



find 



\ait + T)\^\aiT)\^)r^ J2 



^inujMt / in{S<p{t+T)-54i{t)) 



Zo + 2 ^ Z„ cos(nLJMt)e^^^'('''^(*», 



where Zj, ^ \J2r 



-oo "-m^m-n I 



(for more details see 



the supplementary material). This theory explains at 
least qualitatively the shape of the correlator even deep 
in the quantum regime (see Fig. [5]). 

Finally, we can also quantify the decay of g^^\t) in the 
red detuned regime, before the start of the self-induced 
oscillation. Semiclassically, there is no mechanical mo- 
tion in this regime, and thus the dynamics g^'^\i) cannot 
be explained. But we have found that the correlator de- 
cay is then well described by the optomechanical cooling 
rate (see supplementary materials). 

To summarize, in this paper we investigated quantum 
signatures of light and mechanics for an optomechani- 
cal system in the "mechanical lasing" regime. We found 
that, for strong optomechanical coupling [g^ n), for a 
range of detuning and input power the steady state me- 
chanical Wigner density contains strong negative parts, 
signaling stable non-classical states. Back-action evading 
homodyne measurements can be used to reconstruct the 
these Wigner densities. In addition, the two-point pho- 
ton correlator g'^'^\t) displays clear signatures near the 
onset of oscillations. Finally we explained the slow long- 
time decay of the photon correlations as being due to the 
mechanical phase diffusion induced by photon shot noise. 
One should note that experimental observation of some 
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of these photon correlation features does not necessarily 
require being in the nonlinear quantum regime and could 
actually succeed even in present setups. 
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Appendices 



Under the experimentally realistic assumptions, when the 
classical self-induced oscillation starts, it is to a good ap- 
proximation harmonic x{t) ^ x + Acos{ujMt). The laser 
amplitude, influenced by the mechanical oscillation, will 
contain higher harmonics [11] a{t) 



where 



aLJni-goA) 



—nujM + (A + gox) + iK/2 



(A-5) 



and 4>(t) = goAsinliJMt). Here we take the length unit 
to be the mechanical zero point fluctuation xzpf and 
frequency unit to be lom- Jn{x) is the n-th order Bessel 
function. The oscillation amplitude A and equilibrium 
position X can be determined self-consistentl y. W e can 



express '(t) in terms of the solution in Eq. A-5 



7n—-~oo n— — CX3 

oo 

{\a{t + T)\^\a{T)\^}r = Zo + 2j2cos{nujMt)Zn.{A-6) 

71=1 

In this paper we're interested in the strongly quantum 
regime where gM ~ thus in the sideband resolved 



regime we have gu/uJM < 1. From Eq. A-5 we see that 
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A— 0.8n„ . 
decay fit ■ 



The noise spectrum of the radiation force is can be easily 
computed [19] : 



(A- 



The fluctuations of the mechanical oscillator's phase feed 
back to the time-dependence of optical amplitude. Thus, 
under the semi-classical assumption where we take into 
account the photon shot noise but still treat the photon 
amplitude classically, we can modify Eq. |A-6| to be: 



FIG. 6: The decay of the correlation function g''^\t) in the 
red-detuned regime A = —O.Sujm can be qualitatively ac- 
counted for with Fopt obtained from optomechanical cooling. 



only when A ^ xzpf would there be significant contri- 
bution of higher harmonics in the light amplitude a(t), 
which, as can be seen from Eq. |A-6| is also the condition of 
having higher harmonics in g^'^''{t). This explains qual- 
itatively the appearance of higher harmonics when the 
quantum self-induced oscillation gains large amplitude. 

However, even when the self-induced oscillation has 
amplitude much larger than the x^pfi there are impor- 
tant quantum effects that are not accounted for by Eq. 
[5] and Eq. |A-6| As seen in Fig. [Sj the classical solu- 
tion (bottom) is fully periodic, as there is a balance be- 
tween the optical and mechanical dissipation and laser 
driving. However, over the period of 60 cycles, the am- 
plitude of the quantum mechanical g^'^\t) decays signif- 
icantly (top three panels, blue curves). We can qualita- 
tively account for this decay by calculating the influence 
of photon shot noise on the phase fluctuation of the me- 
chanical oscillator. 

Because of the optomechanical coupling in Eq. [T] 
the photon shot noise influences the phase of the me- 
chanical oscillation through radiation pressure F{t) = 
ihujR/L)a{tya{t): 



m) = 

Var((50(t)) = 



(mwM^)^ 4 



dt' F{t') cos((^(t')), (A-7) 

{Sff{,UJm) + SFFi-UJM)) ■ 



\a{t + r)\Mr)\')r= ^ ^« ^ 



n— — 00 



= Zq + 2 Z„ cos{nujMt)e 2 



Var(50(t)) 



Here we assume the phase fluctuation S(j) is Gaussian. 
The result of this semi-classical accounting of the photon 
shot noise can be seen in the red curves in the top three 
panels of Fig. [5] We can see that this simple analysis can 
account qualitatively for the decay of the g^'^\t) in the 
large amplitude self-induced oscillation regime. 

Finally, as we see in Fig. |6j there are also significant 
oscillation structure and decay for the ^'^•'(t) in the so- 
called red detuned regime, before the start of the self- 
induced oscillation. This regime cannot be understood 
at all by the classical picture, since classically the system 
has no dynamics there. The oscillation in t/^^^ (t) can be 
understood as the dynamical response of the mechanical 
oscillator to the quantum fluctuation of the photon field. 
The decay can then be modeled using the theory of the 
optomechanical cooling of the mechanical oscillation in 
the red detuned regime, giving a cooling rate 



opt 



^ZPF 



[Sff{^m) — '5'f_f(— wm)] ■ 



(A-9) 



Here Sf f is the same noise spectrum given by Eq. |A-8| 
As seen in Fig.[6j this rate gives a qualitative account for 
the decay rate of g^'^\t) in the red-detuned regime. 



